PROCEEDINGS 

0F SCIENCE 



Charmless B Decays at CDF 



o 

£ Diego Tonelli* 1 

Fermilab, P.O. Box 500, Batavia, IL, 60510-5011, USA 
E-mail: tonel@fnal.gov 



C3 



The CDF experiment at the Fermilab Tevatron pp collider pioneered the exploration of charmless 
fe-hadron decays in hadron collisions. This provided a unique, rich and highly successful program 
O | of measurements that is currently reaching its full maturity. I report here results on two-body 

decays and the new analysis of B® — > <j><j> decays in an event sample corresponding to 2.9 fb _1 . 
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1. Introduction 

Non-leptonic charmless decays of fe-hadrons are among the most widely studied processes in 
quark-flavor physics. The Collider Detector at Fermilab (CDF) is the first experiment that explores 
these decays in hadron collisions, and has so far unique access to large samples of charmless B° s 
(and A?) decays, necessary supplements to the wealth of precision measurements in the B° and 
B + sectors available from the B-factories [1]. This is made possible by identifying long-lived de- 
cays in the trigger, where dedicated custom electronics [2] reconstructs accurately trajectories of 
charged particles (tracks) in the micro-vertex silicon detector. Tracks in the plane transverse to the 
beam are reconstructed at the first trigger stage (synchronous with the 2.52 MHz bunch-crossing) 
in the multi-wire drift chamber (that covers | tj | < 1) and matched at the second trigger stage (asyn- 
chronous with 25 jUs average latency) with tracks found in the silicon detector. The innermost 
silicon layer, at 2.5 cm radius from the beam, allow measurement of impact parameter (distance 
of closest approach to the beam) with offline-like, 48 /im resolution (including the beam spread). 
This is small with respect to 150 jltm values typical of relativistic ^-hadron decays, allowing dis- 
crimination of long-lived heavy flavor decays from the overwhelming light-quark background. The 
resulting samples of hadronic B decays provided a rich and fruitful program that is currently reach- 
ing its maturity, and whose latest results are shown in what follows. All conjugate modes are 
implied, branching fractions are CP-averages, and K* is shorthand for ^*(892)°. 

2. Updated analysis of B Q S — > 00 decays 

The —^00 decay is mediated by a penguin-dominated b — > sss transition and was first de- 
tected in year 2005 by CDF, with 8 events in a data sample corresponding to a time-integrated 
luminosity of 180 pb _1 [3]. The decay-rate is predicted with large uncertainty and the phe- 
nomenology is enriched by the presence of three different polarization amplitudes, which enhance 
sensitivity to non-standard-model (non-SM) physics. We improved the measurement of branch- 
ing fraction using 2.9 fb _1 of data. The — > 0(— > K + K~)<p(^>- K + K~) rate is measured using 
B® — > J/\j/(^- jU + i u~)0(— > K + K~) decays, collected by the same trigger, as a reference mode as 
follows: 

^(g"^00) jV 00 @{J/y^v, + v,-)e JM 

where N§§ {Nj/^,^} is the number of signal (reference) decays determined in data; the branching 
ratios of the vector mesons are known [4] ; the trigger and selection efficiency to reconstruct a signal 
decay relative to a reference one, is determined from simulation; and the efficiency for 

requiring muon identification on one of the J ' jx\f prongs is accounted for separately by the term, 
determined from inclusive J j\\f in data. We chose B Q S — > J j\\f§ as a reference because it involves 
a decay, and has same number of vertices and final state particles as the signal. The decay 
B° -+ (->■ K + K )K* (-> K+ % ) is even more similar to the signal, but its use as reference carries 
the penalty of the large uncertainty on the ratio of production fractions f s /fd- 

The trigger requires two oppositely-charged particles originated in a point displaced trans- 
versely by the primary pp collision vertex. Trigger tracks are required to have an impact parameter 
of 0.012 < do < 0. 1 cm, to originate at least 200 pLra from the primary vertex in the plane transverse 
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to the beam, and to be azimuthally separated by 2° < < 90°. Different thresholds on transverse 
momenta of particles and their scalar sum (pj > 2 or 2.5 GeV/c, £/?t > 5.5 or 6.5 GeV/c) are 
applied as a function of luminosity to control trigger rates. In the offline analysis, a kinematic fit 
to a common space -point is applied to each track pair. Pairs consistent with being kaons originated 
from <p meson decays are then fit together to a 4-tracks vertex. 

The mass distribution of the four kaons shows a large, featureless flat background. Samples 
of simulated data (to characterize signal, S) and B® mass sidebands (for background, B) are used 
in an unbiased optimization to extract a visible signal. We maximized, over a large set of con- 
figurations of kinematic and vertex quality requirements, the quantity S/^/S + B. This is close to 
optimal for a branching fraction measurement, since S/y/S + B °< l/c%, where a% is the variance 
of a rate in a counting experiment. The most discriminating observables are the transverse decay- 
length of the candidate, Lj(B®) > 330 ^m, and the transverse momentum of the \ower-pj kaon, 
Pt(K) > 0.7 GeV/c. Other used quantities include impact parameter of the highest-px vector me- 
son, c?o(0) > 85 ^m; impact parameter of the candidate, do(B®) < 65 ^m; and the fit-quality of the 
reconstructed secondary vertex, % 2 < 17. 



CDF Run II Preliminary L inl = 2.9 fb'' CDF Run II Preliminary L inl = 2.9 fb" 




m(K*K'K'K") [GeV/c 2 ] m(J/if K'K) [GeV/c 2 ] 



Figure 1: Distribution of K + K K + K mass (left) and J/\j/K + K mass (right) with fit function overlaid. 

The final sample contains 295 ± 20 signal events over a moderate background (fig. 1 , left) dom- 
inated by accidental combinations of random tracks that meet the selection requirements (combina- 
torial). Simulation shows that a 2.5% contamination of B° — > (pK* is also present and contributions 
of the (as yet unobserved) B*' -^K*K* decay are negligible for branching fractions ^(10~ 5 ), as pre- 
dicted by the SM. The events clustering at masses of 5.25-5.30 GeV/c 2 suggest a possible B° — > <p<p 
contribution. A fit where a resonance at the B° mass is allowed to float yields results consistent with 
a background fluctuation of less than 2a, with a corresponding upper limit on £$(B Q —5-00) consis- 
tent with currently known values. This has been included as a systematic effect in the SS{B S — > (pip) 
measurement. The factor 40 gain in signal yield with respect to the previous result [3], in spite of 
just a factor 16 increase in sample size quantifies the effect of the optimized trigger and offline 
selection. Seven layers of silicon sensors extending radially up to 22 cm, and the drift chamber that 
provides 96 measurements between 30 and 140 cm, all immersed in the 1.4 T axial magnetic field, 
provide a mass resolution of approximately 17 MeV/c 2 on the B° — »■ <p<p peak. 

Offline reconstruction of B° s — > J/\j/<p decays is similar to the signal reconstruction, except 
that the 4-track vertex is fit using the known J/y mass as a constraint and at least one track 
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from the candidate J/y is required to have a matching extrapolation outward to track-segments 
reconstructed in the muon detectors (planar drift chambers at |tj | < 1). This requirement selects a 
sufficiently abundant and pure sample rejecting any J/y — > e + e~ contamination. A dedicated opti- 
mization yields 1766 ±48 B® — > J/y<j) events over a low background dominated by combinatorics 
(fig. 1, right) and a residual 4% contribution from mis-reconstructed 5° — > J/yK* decays. Note 
that — > J /yip decays reconstructed from the displaced- track trigger provide an additional 25% 
signal for the analysis of the B^-mixing phase, which suffers from low-statistics issues [5]. 

The relative trigger and selection efficiency is extracted from simulation. Simulated decays 
were reweighed to match the transverse momentum distribution of data because the Monte Carlo 
does not reproduce accurately the observed admixture of triggers with different /^-thresholds. The 
result is Cy/i^/s^ = 0.939 ± 0.009 ±0.030. The first uncertainty is due to finite statistics collected 
by each trigger, the second one arises from the statistical uncertainty in the reweighing. 

The efficiency for identifying a muon is extracted from data using a few thousands inclusive 
J/y —> jU + /x decays as a function of muon pj and in two ranges of pseudorapidity, |tj| < 0.6 
and 0.6 < |tj| < 1.0, corresponding to distinct muon detectors. The J/y decays are collected by 
the same trigger that collects signal decays, which is unbiased for muons, and selected offline as 
follows: Lj(J/y) > 200 /im, % 2 < 5, pj(n) > 1.5 GeV/c, and at least one J/y-tmck that extrap- 
olates outward to a track segment in the muon chambers. The resulting p-p— averaged efficiency is 
£fi = ^/(Nifj. +N2n) = 0.8695 ±0.0044 (stat.), where N[^(N 2 ii) are J/y candidates that have 
exactly one (two) muon(s) identified offline, and the efficiency for identifying di-muons is assumed 
factorizable into the product of single muon efficiencies. 

The dominant systematic uncertainty arises from the unknown polarizations of signal and 
reference decays (6-7%), which impact acceptance. Other contributions include a K/ji difference 
in trigger efficiency due to different ionization in the tracking chamber (3-4%), uncertainty in signal 
mass model (3%), unmodeled backgrounds (3%), /?T-reweighing and background subtraction (1%). 
From the observed number of events, = 295 ± 20 {stat.) ± 12 (syst.) and Nj/^ = 1766 ± 
48 (stat.) ±41 (syst.), we determine 



which, using &(B° S — > J/y$) = (1.45 ±0.46) x 10 3 (the known value updated with the current 
value of f s /f d [4]), yields £g(B° s 00) = [2.40 ±0.21 (stat.) ±0.27 (syst.) ± 0.82(^)] x 10~ 5 . 
The last uncertainty derives from the uncertainty in 33(B° S — > J/y<f>). The result agrees with the 
previous determination, [1.40^q j (stat.) ±0.6 (syst.)] x 10~ 5 [3], and with predictions based on 
QCD factorization, [2.18±0.1^;^] x 10~ 5 or perturbative QCD, [3.53+^ 3 +J ^] x 10~ 5 [6], all of 
which have larger uncertainties. 

This sample of 300 events is already sufficient to explore the rich polarization structure of 
the B® — > <p<p final state. Three states of polarization correspond to the allowed values of orbital 
angular momentum {I = 0, 1, or 2) in the decay of a pseudo-scalar particle into two vector particles 
(B — >• VV). A measurement of their relative contributions provides useful information to the puz- 
zling scenario of B — > VV polarizations. Theory predicts the transverse polarization fraction to be 
small, fj = @{jn\ /trig) 4%, where my (mj) is the mass of the vector (B) meson, but measure- 
ments of fj 50% in B°M _). 0^*(+) decays disfavor this hierarchy [7]. Several ad hoc solutions 



[1.78 ±0.14 (stat.) ±0.20 (syst.)} x 10~ 2 



(2.2) 



8§(B Q s ^J/y<p) 
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may accommodate the discrepancy without invoking non-SM physics. These include enhanced 
annihilation contributions, presence of transversely-polarized gluons, effects of electro-magnetic 
or charming penguins, long-distance re-scattering and others [8]. However, all these explanations 
are either model-dependent or non-conclusive, and the new physics option, e. g., 1 + y 5 terms in 
the amplitude from scalar interactions or super-symmetric particles, remains valid. Polarization of 
the — > is considered crucial to discriminate among models. Ref.[9], for instance, proposes 
to combine it with SU(3) symmetry to constrain the role of penguin annihilation amplitudes. The 
analysis of polarizations is in progress at CDF and 6 (5%) resolutions on amplitudes are expected. 

3. Analysis of #9 s and A° decays into pairs of charmless hadrons 

Two-body decays mediated by the b — > u quark-level transition have amplitudes sensitive to 
7 (or 03), the least known angle of the CKM matrix; significant contributions from 'penguin' 
transitions probe the presence of non-standard-model physics in loops; in addition, the variety of 
open channels with similar final states allow cancellation of many common systematic effects and 
provides information to improve effective models of low-energy QCD. Analysis of these decays 
have a key role in the CDF flavor program, providing unique access to charmless B° s and A° decays, 
and measurements in the B° sector competitive with the B-factories [10] [1 1]. We briefly overview 
here the results based on 1 fb~ 1 of data collected by a dedicated displaced-track trigger similar to 
the one described in the previous section. 

An optimized offline selection isolates a clear signal of a few thousands events over a smooth 
background dominated by combinatorial contributions (fig. 2, left). Absence of event-specific par- 




Invariant rot-mass TGeV/c 2 l K i 

Figure 2: Distribution of 7T7r-mass with fit projections overlaid (left). Separation between K + K~ and K~K + 
states in the space of the dE/dx of the two particles (right). The quantity K is the observed dE/dx around the 
average pion response divided by the difference between average kaon and average pion responses. 

tide identification and limited mass resolution (approximately 22 MeV/c 2 in these decays) cause 
several channels to appear overlapping in a single peak. A five-dimensional likelihood fit rely- 
ing on kinematics differences between decays (correlation between the two-body mass with ar- 
bitrary mass assignment to the final state particles and their momentum imbalance), and particle 
identification from the measurement of specific ionization in the drift chamber (1.5a kaon-pion 
separation from dE/dx, fig. 2, right) allow statistical determination of the individual contribu- 
tions. Large D° — > K~n + samples are used both to model accurately mass line-shapes and for 
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Decay mode Relative branching fraction (33) Absolute 3$(\0 6 ) 

B Q n+ %- %f^~T K + 7C n - ) ) = 0.259±0.017±0.016 5.02±0.33±0.35 

B°^K+K- fxfjfefrf^y = 0.347 ±0.020 ± 0.021 24.4 ±1.4 ±3.5 

B° s ^K-7t + | x = 0.071 ±0.010 ±0.007 5.0 ±0.7 ±0.8 

B°^tc+tc- | x % ( (XTk*k-\ = 0.007 ±0.004 ±0.005 < 1.2 at 90% C.L. 

B°^K+K- = 0.020 ±0.008 ± 0.006 < 0.7 at 90% C.L. 



A° h ^pK- | x S{^r K k') = 0.066 ±0.009 ±0.008 5.6 ±0.8 ±1.5 

A°^p?r- | x = 0.042 ±0.007 ±0.006 3.5 ±0.6 ±0.9 

Table 1: Measured branching fractions. Absolute results were derived by normalizing to the current known 
value 3§(B° -> K+%-) = (19.4±0.6) x 10~ 6 , and assuming f s /f d = 0.276±0.034 and f L /f d = 0.230± 
0.052 for the production fractions [4]. The first quoted uncertainty is statistical, the second is systematic. 



dE/dx calibration. Event fractions for each channel are corrected for trigger and selection effi- 
ciencies (from simulation and data) to extract the decay -rates (tab. 1). We report the first obser- 
vation of the decays B° s — > K~n + , A° — > pK , and A° — > p7T~, and world-leading measurements 
of (upper limits on) the B® — > K + K (B®, — > n + n~) branching fractions [11]. The measured ra- 
tio 3$(t^ b -)• p7T~)/^(A£ -)■ pTT ) = 0.66 ± 0. 14 (stat.) ± 0.08 (syst.) shows that significant pen- 
guin contributions compensate the Cabibbo (and kinematic) suppression expected at tree level. 
We also report first measurements of the direct CP- violating asymmetries Aqp(B q s — > K~7t + ) = 
[39 ± 15 (stat.) ± 8 (syst.)}%, A CP (A° -)■ pw~) = [3 ± 17 (stat.) ± 5 (syst.)]%, and A CP (A° 
pK~) = [37± 17 (stat.) ±3 (syst.)]%. Direct CP violation in the B° — > K + n~ mode is measured as 
Acp(B° ->■ K + n~) = [-8.6 ±2.3 (stat.) ±0.9 (syst.)]%, which is consistent and close to be com- 
petitive with the final B-factories results. The analysis of the 5 fb _1 sample is now in progress. We 
expect to reconstruct more than 15,000 b9% — > h + h'~ decays, where new modes could be observed 
and CP- violating asymmetries measured with doubled precision. 



4. Concluding remarks and outlook 

CDF keeps harvesting from its rich and unique program on charmless B decays in hadron col- 
lisions, pioneered since year 2002. The update of the B® — > 00 analysis to 2.9 fb _1 is the latest 
addition. The 37% uncertainty on the branching fraction is almost half of the previous determi- 
nation, and is dominated by the uncertainty on 3§(B® — > J/\j/(f)), which may soon reduced by the 
Belle experiment. The polarization analysis is in progress with expected uncertainties of 5%. In 
the update to 1 fb _1 of the two-body analysis, the decays B® — > K~7t + , A° — > pn~, and A° — > pK 
were newly observed and their branching fractions and CP- violating asymmetries measured. The 
observed hierarchy 38(A b — > pK~) > 3S(K\ — > pn~) is inconsistent with recent theory predictions 
[12], suggesting significant contributions from penguin amplitudes in hadronic A° decays. 

This fruitful program is reaching its maturity, but it's far from being exhausted. By October 
2010 CDF will have more than 8 fb~ ! of physics quality data on permanent storage, which may 
reach 10 fb~' after one year, if Tevatron Run II will be further extended. These 2.5-10 factors in 
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increase of sample-size over analyses shown here set the scale of what's coming next. Not only 
current analyses, mostly limited by statistical uncertainties and by systematic effects that scale 
with statistics, will be improved. Avenues for new measurements are opening-up: an accurate 
determination of the polarization structure of — > <p<p decays will contribute useful information 
on the puzzling scenario in B — > VV decays; search for a large mixing phase in penguin-dominated 
B® — > (j>(j) decays will supplement the analogous effort in tree-dominated B° — > J/y<f) modes [5]; 
comparison of B° and B® decay rates in common K + n~ final states will provide stringent model- 
independent constraints on non-SM physics contributions [13], exploration of time-dependent CP- 
violating asymmetries in S9 ■> — > h + h ~ decays will increase our knowledge of the angle y (<fo). 

While the first, exciting whimpers of LHC are emerging these days, CDF sits on a goldmine 
of data and a few exciting years of competition with LHCb are coming. 
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